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Abstract. The key to understanding the universal behaviour of systems driven 
away from equilibrium lies in the common description obtained when particular mi- 
croscopic models are reduced to order parameter equations. Universal order param- 
eter equations written for complex matter fields are widely used to describe systems 
as different as Bose-Einstein condensates of ultra cold atomic gases, thermal convec- 
tion, nematic liquid crystals, lasers and other nonlinear systems. Exciton-polariton 
condensates recently realised in semiconductor microcavities are pattern forming sys- 
tems that lie somewhere between equilibrium Bose-Einstein condensates and lasers. 
Because of the imperfect confinement of the photon component, exciton-polaritons 
have a finite lifetime, and have to be continuously re-populated. As photon confine- 
ment improves, the system more closely approximates an equilibrium system. In this 
chapter we review a number of universal equations which describe various regimes 
of the dynamics of exciton-polariton condensates: the Gross-Pitaevskii equation, 
which models weakly interacting equilibrium condensates, the complex Ginsburg- 
Landau equation — the universal equation that describes the behaviour of systems 
in the vicinity of a symmetry-breaking instability, and the complex Swift-Hohenberg 
equation that in comparison with the complex Ginsburg-Landau equation contains 
additional nonlocal terms responsible for spacial mode selection. All these equations 
can be derived asymptotically from a generic laser model given by Maxwell-Bloch 
equations. Such an universal framework allows the unified treatment of various sys- 
tems and continuously cross from one system to another. We discuss the relevance 
of these equations, and their consequences for pattern formation. 



1.1 Introduction 

For a dissipative macroscopic system in thermal equilibrium, relaxation to- 
ward a state of minimum free energy determines the states that the system 
may adopt, and any possible pattern formation. In contrast, if a system is 
driven out of equilibrium by external fluxes, then no such simple description 
is possible, i.e., if a system may exchange particles and energy with multiple 
baths (reservoirs), then the states the system adopts depend not only on the 
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temperatures and chemical potentials of these reservoirs, but also on the rate 
at which particles and energy are injected and lost from the system. This can 
not generally be captured by relaxation to minimise a given energy functional 

Both equilibrium and non-equilibrium systems can be characterised by 
mean-field variables if field fluctuations are negligible (fluctuations can how- 
ever be introduced phenomenologically into the evolution equations if re- 
quired). A mean-field approach leads naturally to the concept of the order 
parameter, and the corresponding order parameter equation. The order pa- 
rameter is either a physical field or an abstract field which acquires a non-zero 
value in the an ordered phase (such as a Bose-condensed or lasing state) , and 
vanishes in the normal state. When considering a spatially inhomogeneous sys- 
tem (with trapping, or inhomogeneous pumping), the order parameter may 
vary in space. When considering a non-equilibrium system, or the dynamics 
of a system as it approaches its equilibrium state, the order parameter may 
also vary in time. In such cases, the order parameter equation describes the 
space and time dependence of the order parameter, accounting for the generic 
features of the system's dynamics. 

One important classification of order parameter equations distinguishes 
whether they describe relaxation towards an equilibrium configuration, or 



phase evolution in a conservative system, or a mixture of the two 32 . For 
a dissipative system, the dynamics may be described by using an energy 
functional written in terms of the order parameter ip and its spatial 

derivatives. The dissipative system dynamics causes this energy functional 
to decay as a function of time, reaching a minimal value at equilibrium, i.e. 
dtip = —Td^J-. The dynamical critical behaviour of such systems has been 



extensively reviewed by Hohenberg and Halperin 14 . Such an approach is ap- 



propriate for many solid-state systems, including in particular non-equilibrium 



superconductivity 21 . In contrast, for isolated systems such as ultracold 
atomic gases, the order parameter obeys conservative dynamics, in which the 
energy functional does not change with time, and the order parameter instead 
follows Hamiltonian dynamics. We will discuss the behaviour of systems that 
lie between these two extremes of purely dissipative and purely conservative 
dynamics, a scenario that includes the non-equilibrium polariton condensate. 

The structure of the energy functional, and thus of the resulting order 
parameter equation, is determined by the symmetries of the order parameter 
space. Taking into account also the fact that near a phase transition, the char- 
acteristic lengthscale of fluctuations diverges, it becomes possible to restrict 
the form of the energy functional by keeping the lowest order derivative terms 
that possess the required symmetries. This makes it possible to divide systems 
into universality classes, depending only on the symmetries and the nature 



of the dynamics 14 16 . Identifying which classes various system belong to 
allows one to draw similarities between systems that are very different in na- 
ture and to predict the behaviour of the new systems that fall into previously 
known universality class. For instance, symmetry under changing the phase 
of of the order parameter restricts the energy functionals to dependence on 
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\tp\ only, and considering the lowest order form that allows for a symmetry 
breaking from disordered (zero) to ordered (nonzero) state gives a potential as 
a quartic polynomial in \ip\. Including the lowest compatible order of spatial 
derivative terms then gives the energy functional whose dissipative dynamics 
correspond to the Ginzburg-Landau equation: 

T= J dVV^-Vr + (p-Uo\M 2 ?, (1-1) 
where \x is the control parameter that forces the system to move from the 



normal state |?/>| = to the ordered state with |i/j| 2 = fi/Uo- Eq. (1.1) is 
expected to be relevant to a physical system in the vicinity of the phase 
transition, where the smallness of the modulus and the derivatives of the 
order parameter allow to keep only the leading order terms in the expansion. 
Higher derivatives and other higher order terms can be kept to allow for more 
complex forms of order and associated phase transitions. 

Understanding the universality class of a given system leads to under- 
standing of fundamentals of the behaviour of that system. The studies of 



vortices and vortex dynamics in superfluid helium 13 led to prediction and 
experimental realisation of vortices first in nonlinear optics [3j , then in atomic 
Bose-Einstein condensates [9] and finally in nonequilibrium solid-state conden- 



sates [26 24 all due the hydrodynamic interpretation of the order parameter 
equations. Spiral waves in biological and chemical systems suggested the ex- 
istence of such meandering waves in class B lasers. Solitary waves in atomic 
systems all have their analogs in nonlinear optics. Finally, much of the exper- 
iments in solid-state condensates are now motivated by finding localised exci- 
tations similar to other system that share the same universality class [T 10 . 
Pattern formation in systems that belong to the same universality class share 
similar properties. Patterns appear in open nonlinear systems when an am- 
plitude distribution of the order parameter becomes unstable above a certain 
threshold. Linear instability gives rise to a so-called pure state that, if non- 
linearities are weak, can dominate the dynamics. Strong nonlinearities may 
mix the eigenvalues leading to various stationary or chaotic combinations of 
pure states with different combinations occupying either all space or different 
space regions. 

In this chapter we shall follow the evolution from equilibrium condensates 
to non-equilibrium condensates to lasers analysing their universality, empha- 
sising similarities and differences. We will discuss in some detail the origin of 
the most general order parameter equation for the laser system, and comment 
on the relation of this order parameter equation to that for cold atoms and for 
non-equilibrium polariton condensation. We will then demonstrate how the 
various terms that may exist in the order parameter equation affect the pat- 
terns which arise, focusing on three cases: the case with homogeneous pumping 
and no trapping, the case with inhomogeneous pumping and no trapping, and 
the case with an harmonic trap. 
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1.1.1 Review of physical systems 

Laser dynamics is described by coupling Maxwell equations with Shrodinger 
equations for N atoms confined in the cavity. Expanding the electric field in 
cavity modes and keeping only the leading order mode leads to the equations 
that couple the amplitude of this mode with the collective variables that de- 
scribe the polarisation and population of the gain medium. Such coupled equa- 
tions are known as Maxwell-Bloch (MB) equations. Lasers are then classified 
depending on the relative order of the loss rates for the electric field, compared 
to the decay rates of the gain medium polarisation and population. The MB 
equations have two homogeneous stationary solutions: nonlasing (zero order 
parameter) and lasing (nonzero order parameter) solutions. The instabilities 
of these solutions, and therefore, pattern formation, are described by univer- 
sal order parameter equations: the complex Swift-Hohenberg (cSH) equation 
for lasers with a fast population inversion and the cSH equation coupled to a 
mean flow if the population inversion is slowly varying. The universal equation 
describing the bifurcation of lasing solution takes the form of a cSH equation 



coupled to a Kuramoto-Sivashinsky equation 23 34 . 

Semiconductor microcavities confine photon modes, which may then inter- 
act with electronic excitations in the semiconductor. If the cavity is resonant 
with the energy to create an exciton (a bound electron hole pair), and if the 
exciton-photon coupling is strong enough, then new normal modes (new quasi- 
particles) arise as hybrids of excitons and photons: polaritons. For low enough 
densities, these may be considered as bosonic quasiparticles, and so can form a 
condensed (coherent) state above a critical density. These are intrinsically non 
equilibrium systems with the steady states set by balance between pumping 
and losses due to the short lifetime of polaritons. Depending on whether the 
emission from the microcavity follows the bare photon or the lower polariton 
dispersion the system shows either regular lasing or polariton condensation 
and in this sense crosses over continuously from weak coupling at higher tem- 
peratures and pumping strengths to strong coupling at lower temperatures 
and lower pumping intensities. Losses in the microcavity systems can be de- 
creased by improving the quality of dielectric Bragg mirrors. The smaller the 
pumping and losses become the closer polariton condensates come to resem- 
ble equilibrium Bose-Einstein condensates (BECs). It seems therefore that 
the unified approach should be possible to describe the transition from nor- 
mal lasers to the equilibrium BECs via polariton condensates. There are some 
other differences between atomic or polariton condensates and normal lasers. 
The operation of a photon laser is based on three ingredients: a resonator for 
the electromagnetic field, an gain medium and a excitation mechanism for 
the gain medium. When excited, the gain medium will undergo stimulated 
emission of radiation that amplifies the electromagnetic field in the cavity. 
In contrast, for polariton condensates there is instead stimulated scattering 
within the set of polariton modes, and condensation can take place without 
any inversion of the gain medium, and thus potentially having a lower thresh- 
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old 15 . A microscopic theory would be required to fully describe how all of 
these aspects cross over from polariton condensation to lasing, however given 
the universality of order parameter equations, one may hope to write a single 
order parameter equation which captures these different regimes by varying 
appropriate parameters. 



1.2 Derivation of order parameter equations 

In this section we show how various order parameter equations arise in descrip- 
tions of lasers, and how these relate to the order parameter equations relevant 
to ultracold atoms and polariton condensates. We start with a mean field 



(semiclassical) model of a laser, the MB system of equations. In section 1.2.2 
we show how the assumption of small relaxation times for atomic polarisation 
in comparison with the cavity relaxation time reduces these equations to the 
complex Ginzburg-Landau (cGL) equations J2[ or the coupled cGL equation 
and the gain medium population dynamics [38|. These models have been ex- 



tensively used to model non-equilibrium condensates 10,44. 18 42 . We discuss 



how the mode selection, in which a particular transverse mode grows fastest, 
is lost in the derivation of these models. Then, in Sec. |1.2.3| we instead follow 
the derivation used in [30] based on the multi-scale expansion technique to 



derive the cSH equations for class A and class C lasers [8]. In section 1.2.4 
we then discuss how nonlinear interactions appear in these equations, and 
discuss the interpretation of these equations as order parameter equations for 
polariton condensates. If the reservoir dynamics is slow in comparison with 
time evolution of the order parameter these equation should be replaced by a 
coupled system explicitly modelling the reservoir dynamics. In the limit of the 
long life-time of the particles the system becomes the Gross-Pitaevskii (GP) 
equation aka the nonlinear Schrodinger (NLS) equation that describes atomic 
BECs. 



1.2.1 Maxwell-Bloch equations for a laser 

We start with the MB equations for a wide-aperture laser with an intracavity 
saturable absorber with multiple transverse modes in the single longitudinal 
mode approximation |27| 
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— - i\7 2 E = P g - P a - (1 + iA e )E, (1.2) 
BP 

T U ^ + (l + iA g )P g = EG, (1.3) 
BP 

T ±a ^ + (l+iA a )P a = EA, (1.4) 

BC 1 

r g -^ = G -G--(E*P g + EP;), (1.5) 

r a —=A -A--(E*P a +EP:), (1.6) 



where the complex field E is the envelope of the electric field, the real functions 
G and A are the population differences for gain and absorption media, the 
complex functions P g and P a are the envelopes of polarisation for gain and 
absorption media. Go and Aq are the stationary values of the population 
difference in the absence of the laser field; they are proportional to the external 
gain and losses in the system. The parameter D = Tj_ a T a [j, 2 / (T± g T g [j, 2 ) is the 
relative saturability of gain and loss media and \x a and fi g stand for the atomic 
dipole momenta. The parameters r± a _ g and r a ^ g are the relaxation times for 
atomic polarisations and population differences scaled by the cavity relaxation 
time, the time t is also scaled by the cavity relaxation time. The parameters 
A a — A e = (ui a — cu c )Tj_ a and A g — A e = (ui g — u! c )Tj_ g are detunings between 
the spectral line centre and the frequency of empty cavity mode u> c . Without 
loss of generality we work in rotating frame such that the the fast time is 
eliminated via introduction of A e . The spatial coordinates are rescaled by the 
width of the effective Fresnel zone. 



1.2.2 Fast reservoir dephasing limit 



Following 35 11 we assume that Tj_ a ,Tj_ g are small and consider the first- 
order approximations to Eqs. ( 1.2||1.6 1. Keeping up to the linear terms in these 
small quantities we get 



Pn 



GE 
1 + iA Q 



T ±9 



(GE) t 



{l + iA g y 



(1.7) 



and similar for P a . The equation on E after we substitute these expressions 
for P g and P a becomes 

(JG 

(1 + i v)g t ~ *(V 2 - A e )e = [(1 - iA g )g - (1 - iA a )a - l]e, (1.8) 



where 



r, = -2r ±g gA g /(l + A 2 g ) + 2r ±a aA a /(l + A 2 a ), 



(1.9) 

and where we rescaled fields as e = E/(l + A 2 g ),g = G/(l + A 2 g ) and a = 
A /(\ + A 2 ). In writing Eq. (1.8) we kept leading order contributions in the 



1 Universality in modelling pattern formation 7 



imaginary coefficient of time derivative (which is of order 0(Tj_ g ,Tj_ a )). The 
real coefficient of the time derivative we kept to 0(1) in rj_ g and rj_ a . The 
equations for scaled gain and absorption media populations to leading order 
take the forms: 

(i + M 2 ) 5 , (1.10) 

(l + d|e| 2 )a, (1.11) 
where a = A /(l + A*), go = G /(l + A 2 g ) and d = D(l + A 2 g )/(1 + A*). 



ag 

T °di =9 ° 
da 

Ta di 



Fast reservoir population relaxation 



The system of equations ( |1.8 1.10[ |1.11| can be simplified under more strin- 
gent restrictions on parameters. In the limit of fast population relaxation times 
r 5 ,r a <l (class A and C lasers) Eqs.( |1.10| |1.11[ ) give 



9 = 



9o 

i + H 



«0 



l + d\e\ 



and Eq. (1.8) becomes 



(l + i77(e))^-zV 2 e = 



(1 - iA g )g a (1 - iA a )a Q 



where the coefficient i] is given by 



77(e) = -2 



TXggAg 



l + d|e| 5 



r± a aA a 



- 1 



e, 



(l + 43)(l + |eP) (l + 4*)(l + d|e|») 



(1.12) 



(1.13) 



(1.14) 



Close to the emission threshold |e| 2 <C 0, which allows a cubic approxima- 
tion for the nonlinear terms we get the complex Ginsburg-Landau equation 
(cGL)|] 



(i - r)(0)) 



de 
di 



-V 2 e + Ve + U\e\ 2 e + i[a - /3|e| 2 ]e, 



(1.15) 



where we let a = g -ao-l, (3 = g -a , U = da a A a ~g A g , V = g A g ~a A a . 

The cGL equation is not a very accurate model of a laser since it does not 
take into account the selection of transverse modes. The lasers emit particular 
transverse modes that depend on the length of the resonator. By making the 
assumption that T± g ,Tj_ a — > 0, we assumed that the gain line is infinitely 
broad. In order to take into account the tunability of lasers that allows spatial- 
frequency selection a more careful derivation of the order parameter equation 
is required, which does not take this limit of fast polarisation relaxation. 
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1.2.3 Multi-scale analysis of the Maxwell-Bloch equations 

In this section we derive the complex Swift-Hohenberg equation capable of 



selecting particular transverse modes from the MB Eqs. (1.2 1.6 1. A similar 
derivation has been done for the MB equations taking into account gain only 
and assuming that A g is small 25 . Here we shall only assume that V 2 — A e 



is small and use it as a small parameter, e(V 2 — A e ). We apply the technique 
of multi-scale expansion to E, P g , a , G and A looking for solutions in the form 
of a power series expansion in e, and introducing two slow time scales T\ = 
et, T 2 — e 2 t, so that d t = £<9ti + e 2 <9r 2 - Next we solve equations at equal 
powers of e. At the leading order we get non-lasing solution (E, P g ,P a , G, A) = 
(0, 0, 0, G , Ac). At 0(e), (E 1 ,P gl ,P al ,G 1 ,A 1 ) = (V>, G a ^/(l+iA g ), A iJj/(1+ 
iA a ), 0,0), where ip is a yet undetermined complex field and Go and Ao are 
linked via 1 = Go/(l + iA g ) — Aq/(1 + iA a ). This condition specifies Go 
and Aq at the threshold for laser emission as G cr it = A a (l + A 2 g )/(A a — A g ) 
and yl cr it = A g (l + A 2 a )/(A a — A g ). We make near-threshold assumption 
G = G crit + e% and A = A CTit + e 2 l a . At 0{e 2 ) we get 

J£ = i{V 2 - zA e )V> + P g2 - P a2 - E 2 , (1.16) 
oil 

BP 

r ±g -^- + (l + iA g )P g2 = E 2 G , (1.17) 
BP 

r± a ^- + (1 + iA a )P a2 = E 2 A , (1.18) 

Q = -G 2 -^p; x +rp 9 i), (1-19) 

Q = -A 2 -^P* al +rPai). (1-20) 
From these equations we get the compatibility condition 

(1 + G Q T± g - A T ±a ) H = i(V 2 - A e )^, (1.21) 

and expressions for P g2 , P a2 , G 2 and A 2 

P 9 2 = -rxgGo^r, P a2 = -r ±a A ^-, (1.22) 

G - Go|Vj|2 A A ° m2 1 (123) 

where we let E 2 = and denoted r± gi ± a = r± g ,± a / (1 + iA 9ia ), Go = Gq/(1 + 
iA g ) and A = A /(l + iA a ). At 0{e 3 ) we get 
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dip 



p, 



93 



\8T 2 
( dP al 
\ dT 2 



Pa3 - 

dP g2 

m. 

dP a2 



E 3 



(1 + iA g )P g3 



(1.24) 

E 3 G +ijj(G 2 + l g ), (1.25) 



+ (1 + iA a )P a3 = E 3 A + iP(A 2 + l a ), 



r g g^ = -G 3 --^p; 2 + rp a 2h 

t„^ = -A 3 - ^-{i>P* a2 + tfPaa). 



9Ti 



(1.26) 
(1.27) 
(1.28) 



The compatibility condition at this order after we substitute ( 1.22 )-( 1.23| 
gives 



dip 



(1 + G a Tj_g - A Tj_ a ) — + fj_g 



o2 ~ dP a2 
- r± a - 



l + iA„ 



1 + iA n 



1 + A g 



dTx 
ApD 
l + Al 



(1-29) 



We use Eqs. (1.21) and (1.22) in Eq. (1.29), collect the derivatives as d t = 
edr 1 +e 2 dT 2 , absorb e into ip and V 2 — A e and replace e 2 l g {e 2 l a ) with Go — G cr i t 
(Aq — A CI it) as expected. The result is the cSH equation 



(1 + G Q Tj_ g - A T± a 



dtp 
~dt 



z(V 2 - A e )iP 



(^l 9 Go 



(i + G T ±g - A T± a y 



;(V 2 



+ jip — 



Go 



A D 
1 + Zi 2 



A G f^ 
(1.30) 



where 7 = (G - G crit )/(1 + iA g ) - (A - A ait )/(1 + iA a ). 

We can simplify the coefficients by considering a limit A g , a <C Tig, a <C 
1, neglecting 0(A 2 a ) and 0(t]_ a A g ^ a ) terms and keeping only the higher 
order terms for real and imaginary parts of the coefficients. This leads to the 
following general form of the cSH equation 



dip 
~dt 



i(V 2 - A e )i> - S(V 2 - A e ) 2 ip 
+ (a-iV)ip-(p + iU)\ip\ 2 ip, 



(1.31) 



with the energy relaxation r\ = —2GoA g Tj_ g + 2AoA a Tj_ a , the coefficient of 
superdiffusion 5 — r]_ g Go — r^_ a A a , the effective pumping a = Go — Aq — 1, 
the effective repulsive potential V — GqA 9 — AoA a , the cubic damping f3 = 
Gq — A D and interaction potential U — A DA a — GoA g . 

Apart from nonlinear optics and lasers the cSH equation provides a re- 
duced description of a variety of other systems [8] , such as Rayleigh-Bernard 
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convection [7], Couette flow 28 , nematic liquid crystal [5], magneto convec- 
tion [7] and propagating flame front 29 among others. 

Similar to other universal equations the cSH equation can be derived phe- 
nomenologically from general symmetry considerations. Assuming that the 
system is characterised by an instability at k c 7^ 0, the dominating growth 
rate (Lyapunov exponent) can be approximated close to k c by a parabola 
that takes positive values in the neighbourhood of k c . To the lowest degree of 
approximation this can be modelled by 

X = a-S(k 2 -k 2 c f+i(k 2 -k 2 c ), (1.32) 

where a is a control parameter that takes Re(X) into the positive range of 
values. A linear model that has the corresponding dispersion has to be com- 
plemented with a nonlinear term in order to prevent the infinite growth of 
unstable modes. The simplest form of such nonlinearity that preserves the 
invariance of the field phase is the cubic nonlinearity |V'| 2 '0- So the minimum 
equation that describes a class of phenomena in nonlinear optics in the lowest 



order approximation coincides with the cSH equation (1.311 



Slow population evolution 

For wide aperture CO2 and semiconductor lasers the cSH equation introduced 
in the previous section is not a good model. The population dynamics is slow 
which corresponds to the case of the stiff MB equations that occurs when 
the parameter b = r_L S)0 /r gja , that measures the ratio of the polarisation 
dcphasing to the population deenergisation rate, becomes small. The order 
parameter equation in this case is not a single equation and the analysis of 



the previous section should be revised taking into account smallness of b 25 
Instead of going through the multi-scale analysis we note that we can consider 
gain selection separate from population evolution and therefore rewrite Eqs. 
( 1.8|1.10|l.ll| to include the gain selection mechanism 



(1 + iri) 9 ^ = *(V 2 - A e )i, - S(V 2 - A e ) 2 ^ - ^ 
at 



'a 

lg — -G -(1 + M 2 )G, (1.34) 



[(l-iA g )G-(l-iA a )A]1>, (1.33) 
dG 
' dt 

r) A 

r a C ^=A -(l + D\^\ 2 )A (1.35) 
One may note that in the limit that r„, r a are small, this equation reduces to 



Eq. (1.31) 



1.2.4 Modelling exciton-polariton condensates 



The cSH order parameter equation derived above from the MB equations of 
a laser can also describe the polariton condensate. In this section we discuss 
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how such an equation can arise for the polariton system, and the meaning 
the various terms would acquire in this context. We also make contact with 
the limiting cases which correspond to ultracold atomic gases. For the polari- 
ton condensate we interpret ip in Eqs. ( 1.33^ 1.35 ) as a scalar mean field of 
a polariton matter-wave field operator <P(r,t). We begin by considering the 
basic energy functional for a polariton condensate. In addition to the kinetic 
energy, and any external trapping potential, one must also take into account 
repulsive interactions of polaritons. These interactions predominately come 
from the short ranged electron-electron exchange interactions (when two ex- 
citons swap their electrons). This interaction gives rise to a cubic nonlinear 
term — iUo\ip\ 2 tp just as in the right-hand side of Eq. (1.31 ). Rather than cou- 



pling the order parameter equation to the dynamics of the gain medium, one 
should instead consider coupling of the order parameter equation to the equa- 
tion describing the density of noncondensed polaritons (reservoir excitons), 

G ~ 



20 42 



that may also contain a diffusion term. 
In the limit of vanishing gain and losses, the order parameter equation 
becomes the NLS equation also used to model a Bose-Einstein condensation 
of ultracold atoms: 



dtp 



iV 2 4> + iV{r)i> - iU \ip\ 2 ip. 



(1.36) 



For an ultracold atomic gas this equation can also be derived microscopically 
from the Heisenberg representation of the many-body Hamiltonian using the 
language of second quantisation. For the case of an ultracold atomic gas, one 
may also include the dissipation that arises from collisions of condensate atoms 
with non condensed thermal cloud in this equation. This process leads to en- 
ergy relaxation and atom transfer between the condensate and the thermal 
cloud. This can be modelled by writing the quantum Boltzmann equation, i.e. 



kinetic equation, describing the dynamics of the populations of states 31 12 



The net rate of atom transfer rj as the result of such collision can be repre- 



sented by replacing the time derivative in Eqs. (1.36) as d t — > (1 + irj)dt. This 
parameter 77 depends on the temperature and the density of the noncondensed 
cloud. Similar mechanism of energy relaxation exists in polariton condensates 
and have been phenomenologically introduced into various models of polari- 



ton condensates 44 43 . Note that such energy relaxation follows directly 



from the MB equations as indicated by Eqs. (1.8 1.31 1. The interactions with 



noncondensed cloud may enhance this coefficient. 

In addition to the terms mentioned so far, the polariton system differs 
from ultracold atoms, but is similar to the laser, in having also terms describ- 
ing gain and loss, i.e. pumping and decay. Including these terms, and allowing 
them to potentially depend on wavevector, gives a modified cSH model that 
includes all possible previously discussed limits of lasers, nonequilibrium po- 
lariton condensates and equilibrium atomic BECs: 
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(l+«7) 



dip 
0t 



z(V 2 - A e )4> - <5(V 2 - A e fi, - V 



+ [(1 - zZ\ 9 )G - (1 - i^ )]V» - *wiv, 
r „^ =Go -(l 



9 at 

' at" 



1 2 )G, 
A -(l + £>|Vf)A 



(1.37) 
(1.38) 

(1.39) 



Note that = GA g — AA a gives rise to a reservoir potential which causes 
the blue-shift in the condensate 



39 



Some limiting cases of Eqs. (1.37-1.39) have been previously considered. 
Assuming S — > and fast relaxation of reservoirs (t , t 9 — > 0) leads to the cGL 
equation introduced for polariton condensates in [18] . In the limit 5 — > and 
assuming the slow relaxation of the noncondensed reservoir gives rise to the 



model of atom laser 20 that has proved effective for polariton condensates 



42 22 33 44 . Finally, in the limit of vanishing losses and gain all systems 



approach the conservative NLS equation. 



1.3 Pattern formation and stability 

Having discussed the physical origin of the order parameter equations of po- 
lariton condensates, lasers and atomic condensates, this section now discusses 
the consequences of the form of the order parameter equations for pattern 
forming and stability analysis. We discuss three cases: the entirely homoge- 
neous case, the case in which the pumping (gain) is localised, and the case in 
which there is inhomogeneity of the condensate mode energy (i.e. trapping) 
as well as pumping. The homogeneous case is most relevant to wide aperture 
lasers with electrical pumping. For polariton condensates and photon conden- 
sates with external pumping, the second and third cases are more relevant. As 
one goes toward equilibrium systems (such as atomic condensates), the role of 
trapping potentials to confine the condensate becomes more important, and 
so the third case is most relevant in this limit. 

All three cases can be written as the short population relaxation time limit 



of Eq. (1.371 



(1 +ir>)~ = («(r) - m 2 )^ + *(V 2 - V(r) - U \^ 

+2SA e V 2 iJj - <5VV, (1-40) 
but we will rescale the equation in different ways for the different cases. 

1.3.1 Behaviour of homogeneous order parameter equation 

We begin by reviewing the simplest case, of linear stability analysis about 



the uniform solution p — ip 2 — a/(/3 + Uoij) in Eq. (1.40). This uniform 
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solution should be a stable solution as long as SA e < 0. One may then consider 
perturbations of the Bogoliubov-de Gennes form: 

i/j = U + ue l(k ' x - wt) + v * e -i(*-*-"*t)\ e -*M* (i.4i) 



where the chemical potential is fx = Uop. For this ansatz to solve Eq. (1.40) 
(at linear order order in u,v) requires that: 

p 2 (U% + /3 2 ) = |S| 2 + w(l - iv)S* - u(l + iv)S - uj 2 {1 + v 2 ), (1.42) 

where S = (fc 2 + pUo) + i(2SA e k 2 + 8k A + pj3) (making use of the steady state 
values of p and p). For an equilibrium condensate a = j3 = S = r) = one can 
recover the expected Bogoliubov spectrum from Eq. ( |1.42 1 : 



w B (fc) = Vfc 2 (2pC/o + fc 2 )- (1-43) 

Alternatively, in the cGL regime (5 = 0) with r\ — one recovers the dissipa- 
tive spectrum obtained previously (4l) : 

w C GL(fc) = -ia ± V^s(fc) 2 - a 2 , (1.44) 

which is imaginary for small k 7 and then above a critical k (set by ws(fc) = a), 
the imaginary part becomes a constant —a and a real part appears. Introduc- 
ing the remaining terms gives 

UcSH(k) = ——3 [-i(a + Xfc) ± V-" 2 + tk(2p(Uo - W) + ejOl ( L45 ) 

1 + T] Z l J 

where = fc 2 [?7 + <5(2Zi e + fc 2 )] and e fc = fc 2 [l - 7]5(2zi e + fc 2 )]. Note that for 
fc = 0, one always has a mode at zero frequency, as expected given the phase 
symmetry breaking present in the ordered phase. As long as 7/ > 0, 5A e > 0, 
the imaginary part grows for large fc, since such a case describes pumping 
that suppresses high energy (momentum) modes. If rj + 25 A e > 0, the modes 
are always decaying, but if A e < —rj/25, it becomes possible for the CSH 
term to make the uniform part unstable — the exact critical A e depends in 
a non-trivial way on the remaining parameters. Other instabilities may also 
arise due to the content of the square root term. 



1.3.2 Inhomogeneous pumping 

We next consider the effect of inhomogeneous pumping, comparing the be- 
haviour of cSH equation and cGL equations when used to model polariton 
condensates. As the first example we consider a small pumping spot. This 
geometry has been considered extensively in experiments (40,39,6,33] and 



theory 44 . Our starting point is to consider Eq. (1.401 that we rewrite as 
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(l + ir,(P))^ = ( P(v) - -.. - A/'d-lk- ) r + ,(V- - - HP) - I , t\), ■ 



-28A e V 2 ip-6V 4 i>,. (1.46) 



We take 



P(r) = 4exp(-0.05r 2 ), 7c = 0.3, A = 0.075, 

7]{P) = 0.025P(r), V(P) = 1.25exp(-0.45r 2 )P(r). (1.47) 

In writing the last two expressions we recalled that rj, the energy relaxation 
parameter representing the rate of transfer between the noncondensed and 
condensed polaritons, depends on the density of the noncondensed cloud. We 
also assumed spatially dependent energy shifts coming from strong mutual 
repulsion |17| , so that the repulsive force coming from potential V varies with 
density of the condensate. 

We compare two cases: the cGL equation by letting 5 = in Eq. ( 1.46| and 



the cSH equation with S = 0.1, A e = —1. In the case of the evolution accord- 
ing to the cGL equation the system reaches the steady state, see Fig |l.l[ a) 
which shows on the tomography image Fig |l.l[ c) as a single energy level. The 
evolution according to the cSH equation leads to periodic oscillations of the 
density profile shown on Fig. |l.l[ b). The corresponding tomography image 
on Fig. 1 1 , 1 1 shows several discrete energy levels. Similar behaviour has been 
observed in some experiments, eg. [33] . 



1.3.3 Inhomogeneous energy (trapping) 

We now turn to consider the behaviour in the presence of an harmonic trap [4j 
19 . We will consider how the presence of the dissipative terms in the general 
order parameter equation affects the stability of known solutions of the Gross 
Pitaevskii equation. We will look both at linear stability analysis (where one 
can gain insights from analytical results found by considering the perturbative 
effect of dissipation), as well as full numerical time evolution to find the new 
steady states. 

As a starting point in the absence of dissipative terms, the Gross-Pitaevskii 
equation: 

d t i) = i(V 2 - r 2 - L/oM 2 )V (1-48) 

can be approximately solved by the stationary Thomas-Fermi profile with 
dtip = —ifiip and \ip\ 2 = 0(fi — r 2 )(/i — r 2 )/U~o. This density profile results 
from neglecting the kinetic energy. This is valid as long as the cloud size 
r TF = is large compared to the healing length 1/y/Uo, i.e. for y/Uo/j, ~> 1. 
This stationary Thomas-Fermi profile gives a simple prescription for how to 
find the density profile in a given potential landscape. However, as will be 
discussed below, the stationary profile does not necessarily remain stable in 



the presence of the additional terms in Eq. (1.401 
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Fig. 1.1. Comparison between solutions of the cGL and cSH equations for a single 
pumping spot. Time evolution of the density measured across the pumping spot 
(a) for the cGL equation and (b) cSH equation. Energy of the solutions across the 
pumping spot (c) for the cGL equation and (d) cSH equation. The results of the 
numerical simulations of Eq. (1.46 1 with the parameters given in (1.471. 



Stability analysis 



Starting from Eq. (1.40) with a, (3,6, r\ — 0, we consider in turn the effects 
introduced by adding these dissipative terms. We restrict to considering A e < 
and neglect the superdiffusion term; after rescaling parameters, we may 
write: 



/3|V>| 2 + S(2A e ~ V 2 )V 2 - 2ifjd t 



1> 
(1.49) 

in which all dissipative terms are placed on the right hand side. We then 
proceed by considering normal modes around the stationary solution in an 
approximation where the quantum pressure terms can be neglected, this is 
done by writing the equations in terms of density and phase and neglecting 



16 



N.G Berloff and J. Keeling 



all quantum-pressure-like terms: 

dtp + v • (pV0) = (& - ~p P + 2fjd t <f> - 25(V0) 2 ) P , (1.50) 

2d t ct> + (V0) 2 + r 2 + p = 5(2A e V 2 ~ V 4 )0. (1.51) 

In the absence of the dissipative terms, this problem is the two-dimensional 

Linearising these equations us- 

he~ 



analog of that studied by Stringari 36 37 
ing p — > p + he~ lut ,4> — > <fi + ipe~ %u)t yields normal modes with frequen- 



cies u) r . 



2(s + l)n + s and density profiles given by hypergeomet- 



ric functions h(r, 9) oc %Fx{— n, n + s + 1; s + 1, r )e 



is6 „s . 



here n is a radial 



quantum number, and s is an angular quantum number. Including dissipative 
terms, these normal mode frequencies acquire imaginary parts, describing ei- 
ther growth or decay of such fluctuations. Instability of the stationary state 
occurs when at least one of these fluctuation modes grows. 

To account for dissipative terms perturbatively, it is enough to take the 
normal mode functions found in the absence of dissipation and find the first 
order frequency shift induced by the dissipative terms. At first order in the 
dissipative terms, there is no change to the density profile; however a non- 
zero phase gradient V(/> does appear at first order in the dissipative terms. 
Vanishing of the current at the edge of the cloud then requires p = 3<S/(2/3 + 
3fj). 



Formally one may write the linearised form of Eq. (1.50) in the form 
—iiOs.nXs.n( r ,Q) = (C 1 * * 1 + C^ 1 ') Xs.n m which C is a matrix of differential 
operators acting on the fluctuation term \ = (h, v) T - By identifying the dissi- 
pative terms as CP-\ standard first order perturbation theorjj^] then yields the 

first order correction: uifS = i{Xn) , CP^Xn)) /(Xn°i , Xns) where angle brackets 
indicate the appropriate inner product. Following this procedure, one eventu- 
ally finds 



,(i) - 



i 

2N 



2irrdr 



(h^f (c 



■qp 



(20 + ti)n) +5h® 



4 - ^V 2 
2 



v 2 /4°2 

(1-52) 

where the normalisation N — j 2irrdrh 2 ls and integration is over the area 
of the Thomas-Fermi profile r 2 < p . The hypergeometric form of the zero 



order functions hns allows Eq. (jl -52h to be evaluated analytically. The terms 



proportional to S in fact vanishes, and the remaining term can be written 
(making use of the above value of p) as: 



4/3 + 6^ 



(6/3 + 3fj) 



s 2 + (c^) 2 



•2(w<2?) 2 , 



4/3 -3?7 



(1.53) 



3 Some care must be taken since the operator C is not self adjoint and so the 
left and right eigenstates of must be found separately. This is easiest if one 
replaces the variable ip by u r = d r Lp,ug — (l/r)detp = (is/r)ip, in this case the 
right eigenstates (h,u r ,ue) corresponds to the right eigenstate (h,2pu r ,2pug). 
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Crucially, (w„ s ) 2 as given above grows only linearly with s. Thus, at large s 
the ratio in parentheses tends to one, and so u>„s — ► W&/(2j3 + 3fj) > This 
positive value means that there is an instability, even for non-zero fj. Although 
neither fj nor 5 remove the instability in this perturbative approach, this does 
not prevent these terms from restoring stability via higher order corrections. 
This needs to be checked by numerical simulations. 



Vortex lattices 



Having seen that the stationary profile possesses an instability, we next con- 
sider the behaviour resulting from this instability. In order to reach a final 
configuration, it is necessary to restrict the pumping to a finite range. We 
thus take and <5(r) = a 8(r — r), So = 8, /3 — 0.6, where the pumping is 
a flat top of the radius ro = 7. For fj = 6 = 0, this model has been found 
to evolve to a rotating vortex lattice [18]. If one also includes the superdiffu- 
sion present in the cSH model, one finds that (in contrast to the linearised 
analysis) this may arrest the instability to vortex formation, and thus lead 
instead to an oscillating vortex-free state. Fig . 1 1 . 2 1 compares the profiles that 
result from the numerical simulation of Eq. ( |1.46[ ) for the cases of the cGL 
equations for fj = (Fig. |1.2[ a)), fj = 2 ( Fig. |l.2| (b)) and the cSH equation 



with fj = 5 = 0.2 and A e = -0.5 (Fig. 1.2 (c)) 



Although the presence of fj does not remove the instability, it does sig- 
nificantly effect the resulting rotating profile. This can both be seen in the 
numerical results show n in F ig. |1.2[ and can also be understood by considering 
the (5 = limit of Eq. ( 1.50 ), written in a rotating frame. In a rotating frame, 



we consider solutions to Eq. (1.49) of the form 



(p - 2iil r xV)ip 



(1.54) 



such that the time dependence has two parts: rotation with angular velocity 
$1 = (0,0, J?), and phase accumulation at rate p. Inserting this ansatz into 
Eq. ( 1.49 ) and then making the Madellung transform, with neglect of quantum 



pressure terms, one finds: 

V • [p(V<f> - n x r)] = a- f3p-fj(n + 2Q x r ■ V0) p (1.55) 
-p + (V0 - n x r) 2 + (1 - f2 2 )r 2 + p = (1.56) 

These equations can be satisfied by setting V0 ~ 1~2 x r which yields: 

= a - j3p - fj (p + 2Q 2 r 2 ) , -p + (1 - f2 2 )r 2 + p = (1.57) 

These give two equations for p which are both satisfied if: 



p = p-(l- n 2 )r 2 = j [a - fj(p + 2f2 2 r 2 )] 



(1.58) 
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Fig. 1.2. Density plots of the polariton condensate in a harmonic trap according to 
the cGL equation with 77 = (a), r\ — 0.2 (b) and the cSH equation with S — r\ — 0.2 
and A e = —0.5 (c). Top panels: luminosity is proportional to density. Red circles 
show the size of the pumping spot. Bottom panels: (solid lines) density of the cross 
section of the condensate at various times; (dashed lines) analytic approximations, 
given for (a) and (b) by Eq. (1.581, and for (c) by the Thomas-Fermi profile 0(n — 



r 2 )(n - r 1 ) with fi = 3a /(2/3 + 3fj). 



hence fi 2 = j3/(fi + 2rj),li = 01 /(fi + v)- This indicates that while for 77 = 0, 
the lattice rotates at Q = 1, cancelling out the trapping potential, for finite 
77, the rotation velocity decreases, hence the density profile becomes non-flat, 
as seen in Fig. |1.2| 

In the above, V</> ~ O x r would require the phase profile to mimic solid 
body rotation. For a condensate, this cannot be exactly satisfied, but can 
be approximately satisfied (on a coarse grained scale) by having a density of 
vortices fl. Since increasing f] causes fi to decrease, a sufficiently large value 
of fi can in effect kill any finite vortex lattice by reducing the vortex density 
to values so that the number of vortices falls below one. 



1.4 Conclusions 

We reviewed the connection between lasers, polariton condensates and equilib- 
rium Bose condensates from a common framework based on order parameter 
equations. The cSH equations derived for lasers should be applicable to polari- 
ton condensates in the limit of non-negligible interactions and the stimulated 
scattering between polariton modes. The pattern formation in the framework 
of the cSH equations have been well-studied for lasers indicating a wealth 
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of dynamics and phenomena. Some of these phenomena may be achieved in 
polariton condensates. At the same time the stronger nonlincaritics and dif- 
ferent external potentials (engineered or due to disorder) may lead to novel 
properties of the system exhibiting effects not seen in normal lasers. 
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